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Abstract. We discuss notions of Gauss curvature and mean curvature for polyhedral 
surfaces. The discretizations are guided by the principle of preserving integral relations 
for curvatures, like the Gauss/Bonnet theorem and the mean-curvature force balance 
equation. 
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The curvatures of a smooth surface are local measures of its shape. Here we con- 
sider analogous quantities for discrete surfaces, meaning triangulated polyhedral surfaces. 
Often the most useful analogs are those which preserve integral relations for curvature, 
like the Gauss/Bonnet theorem or the force balance equation for mean curvature. For sim- 
plicity, we usually restrict our attention to surfaces in euclidean three-space E'^, although 
some of the results generalize to other ambient manifolds of arbitrary dimension. 

This article is intended as background for some of the related contributions to this 
volume. Much of the material here is not new; some is even quite old. Although some 
references are given, no attempt has been made to give a comprehensive bibliography or 
a full picture of the historical development of the ideas. 



1. Smooth curves, framings and integral curvature relations 

A companion article ISulOSII in this volume investigates curves of finite total curvature. 
This class includes both smooth and polygonal curves, and allows a unified treatment of 
curvature. Here we briefly review the theory of smooth curves from the point of view we 
will later adopt for surfaces. 

The curvatures of a smooth curve 7 (which we usually assume is parametrized by 
its arclength s) are the local properties of its shape, invariant under euclidean motions. 
The only first-order information is the tangent line; since all lines in space are equiva- 
lent, there are no first-order invariants. Second-order information (again, independent of 
parametrization) is given by the osculating circle; the one corresponding invariant is its 
curvature k = 1/r. 
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(For a plane curve given as a graph y = f{x) let us contrast the notions of cur- 
vature K and second derivative /". At a point p on the curve, we can find either one by 
translating p to the origin, transforming so the curve is horizontal there, and then look- 
ing at the second-order behavior. The difference is that for curvature, the transformation 
is a euclidean rotation, while for second derivative, it is a shear {x, y) {x,y — ax). 
A parabola has constant second derivative /" because it looks the same at any two points 
after a shear A circle, on the other hand, has constant curvature because it looks the same 
at any two points after a rotation.) 

A plane curve is completely determined (up to rigid motion) by its (signed) curva- 
ture k{s) as a function of arclength s. For a space curve, however, we need to look at the 
third-order invariants; these are the torsion r and the derivative k', but the latter of course 
gives no new information. Curvature and torsion now form a complete set of invariants: a 
space curve is determined by k{s) and t(s). 

Generically speaking, while second-order curvatures usually suffice to determine a 
hypersurface (of codimension 1), higher-order invariants are needed for higher codimen- 
sion. For curves in E"^, for instance, we need d — 1 generalized curvatures, of order up 
to d, to characterize the shape. 

Let us examine the case of space curves 7 C E'' in more detail. At every point 
p G 7 we have a splitting of the tangent space TpE'^ into the tangent line Tpj and the 
normal plane. A framing along 7 is a smooth choice of a unit normal vector A^i, which 
is then completed to the oriented orthonormal frame (T, Ni,N2) for E'^, where N2 = 
T X Ni. Taking the derivative with respect to arclength, we get a skew-symmetric matrix 
(an infinitesimal rotation) that describes how the frame changes: 




Here, T'{s) = ^ KiNi is the curvature vector of 7, whUe r measures the twisting of the 
chosen orthonormal frame. 

If we fix iVi at some basepoint along 7, then one natural framing is the parallel 
frame or Bishop frame IIBis75ll defined by the condition t = 0. Equivalently, the vec- 
tors Ni are parallel-transported along 7 from the basepoint, using the Riemannian con- 
nection on the normal bundle induced by the immersion in E'^ . One should note that this is 
not necessarily a closed framing along a closed loop 7; when we return to the basepoint, 
the vector Ni has been rotated through an angle called the writhe of 7. 

Other framings are also often useful. For instance, if 7 lies on a surface M with unit 
normal v, it is natural to choose A^i = v. Then N2 = rj T x u is called the cornormal 
vector, and [T^v^rf) is the Darboux frame (adapted to 7 C M C E"^). The curvature 
vector of 7 decomposes into parts tangent and normal to M as T' = k„i^ + Hgrj. Here, 
K„ = —v'-T measures the normal curvature of M in the direction T, and is independent 
of 7, while Kg, the geodesic curvature of 7 in AI, is an intrinsic notion, showing how 7 
sits in M, and is unchanged if we isometrically deform the immersion of M into space. 
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When the curvature vector of 7 never vanishes, we can write it as T' = kN , 
where is a unit vector, the principal normal, and k > 0. This yields the orthonormal 
Frenet frame (T, B), whose twisting t is the torsion of 7. 

The total curvature of a smooth curve is Jk ds. In ISulOSII we review a number of 
standard results: For closed curves, the total curvature is at least 27r (Fenchel) and for 
knotted space curves the total curvature is at least 4tt (Fary/Milnor). For plane curves, 
we can consider instead the signed curvature, and find that j nds is always an integral 
multiple of 2tt. Suppose (following Milnor) we define the total curvature of a polygonal 
curve simply to be the sum of the turning angles at the vertices. Then, as we explain 
in ISulOSI . all these theorems on total curvature remain true. Our goal, when defining 
curvatures for polyhedral surfaces, will be to ensure that similar integral relations remain 
true. 



2. Curvatures of smooth surfaces 

Given a (two-dimensional, oriented) surface M smoothly immersed in E'^, we understand 
its local shape by looking at the Gauss map : M — > given by the unit normal vector 
V = Vp ?il each point p e AI . The derivative of the Gauss map at p is a linear map 
from TpM to T^J^^- Since these spaces are naturally identified, being parallel planes 
in E'^, we can view the derivative as an endomorphism —Sp : TpM TpM. The map Sp 
is called the shape operator (or Weingarten map). 

The shape operator is the complete second-order invariant (or curvature) which de- 
termines the original surface M. (This statement has been left intentionally a bit vague, 
since without a standard parametrization like arclength, it is not quite clear how one 
should specify such an operator along an unknown surface.) Usually, however, it is more 
convenient to work not with the operator Sp but instead with scalar quantities. Its eigen- 
values Ki and K2 are called the principal curvatures, and (since they cannot be globally 
distinguished) it is their symmetric functions which have the most geometric meaning. 

We define the Gauss curvature K := K1K2 as the determinant of Sp and the mean 
curvature H := ki + K2 as its trace. Note that the sign of H depends on the choice 
of unit normal i/, so often it is more natural to work with the vector mean curvature (or 
mean curvature vector) H := Hi'. Furthermore, some authors use the opposite sign on Sp 
and thus H, and many use H = {ki + K2)/2, justifying the name "mean" curvature. Our 
conventions mean that the mean curvature vector for a convex surface points inwards (like 
the curvature vector for a circle). For a unit sphere oriented with inward normal, the Gauss 
map ly is the antipodal map, Sp — I, and H = 2. 

The Gauss curvature is an intrinsic notion, depending only on the pullback metric on 
the surface M, and not on the immersion into space. That is, K is unchanged by bending 
the surface without stretching it. For instance, a developable surface like a cylinder or cone 
has K = because it is obtained by bending a flat plane. One intrinsic characterization 
of K{p) is obtained by comparing the circumferences of (intrinsic) e-balls around p 
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to the value 2tt£ in E^. We get 



= 1 - '-K{p) + 0(£3). 



2776 



Mean curvature, on the other hand, is certainly not intrinsic, but it has a nice varia- 
tional interpretation. Consider a variation vectorfield ^ on M; for simplicity assume ^ is 
compactly supported away from any boundary. Then H = —5 Area /(5 Vol in the sense 
that 



With respect to the inner product r;) :— J £,p ■ f]p dA on vectorfields, the vector 
mean curvature is thus the negative gradient of the area functional, often called the first 
variation of area: H = — V Area. (Similarly, the negative gradient of length for a curve is 
its curvature vector kN.) 

Just as K is the geometric version of second derivative for curves, mean curvature is 
the geometric version of the Laplacian A. Indeed, if a surface M is written locally near a 
point p as the graph of a height function / over its tangent plane TpM, then H{p) = A/. 
Alternatively, we can write H ~ Vm ■ ~ Aj\/x, where x is the position vector in E'^ 
and A m is the Laplace/Beltrami operator, the intrinsic surface Laplacian. 

We can flow a curve or surface to reduce its length or area, by following the gradient 
vectorfield kN or Hiy; the resulting parabolic (heat) flow is slightly nonlinear in a natural 
geometric way. This so-called mean-curvature flow has been extensively studied as a geo- 
metric smoothing flow. (See, among many others, IIGH86llGra87ll for the curve -shortening 
flow and BBraTSI lHui84l |Iirn94i ,Eck04. , WhiOl for higher dimensions.) 

3. Integral curvature relations for surfaces 

For surfaces, we will consider various integral curvature relations that relate area integrals 
over a region D C A/ to arclength integrals over the boundary 7 dD. First, the 
Gauss/Bonnet theorem says that, when D is a disk. 



Here, rfx = T ds is the vector line element along 7, and 77 = T x is again the cornormal. 
In particular, this theorem implies that the total Gauss curvature of D depends only on a 
collar neighborhood of 7: if we make any modification to D supported away from the 
boundary, the total curvature is unchanged (as long as D remains topologically a disk). 
We will extend the notion of (total) Gauss curvature from smooth surfaces to more general 
surfaces (in particular polyhedral surfaces) by requiring that this property remain true. 

Our other integral relations are all proved by Stokes's Theorem, and thus require 
only that 7 be the boundary of 13 in a homological sense; for these D need not be a disk. 
First consider the vector area 



6^ Vol 




6^ Area = 
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The right-hand side represents the total vector area of any surface spanning 7, and the 
relation shows this to depend only on 7 (and this time not even on a collar neighborhood). 
The integrand on the left-hand side depends on a choice of origin for the coordinates, but 
because we integrate over a closed loop, the integral is independent of this choice. Both 
sides of this vector area formula can be interpreted directly for a polyhedral surface, and 
the equation remains true in that case. We note also that this vector area is one of the 
quantities preserved when 7 evolves under the Hasimoto or smoke-ring flow 7 — kB. 
(Compare lLP94llPSW07llHof08l .) 

A simple application of the fundamental theorem of calculus to the tangent vector 
of a curve 7 from pto q shows that 

rq ri 
T{q)~T{p)= / T'{s)ds^ / KNds. 
J p J p 

This can be viewed as a balance between elastic tension forces trying to shrink the curve, 
and sideways forces holding it in place. It is the key step in verifying that the vector 
curvature kN is the first variation of length. 

The analog for a surface patch D is the mean-curvature force balance equation 



T]ds = -fvxdyi^ / / HvdA^ II YidA. 

7 J7 J JD JJd 

Again this represents a balance between surface tension forces acting in the conormal di- 
rection along the boundary of D and what can be considered as pressure forces (espcially 
in the case of constant H) acting normally across D. We will use this equation to develop 
our analog of mean curvature for discrete surfaces. 

The force balance equation can be seen to arise from the translational invariance of 
curvatures. It has been important for studying surfaces of constant mean curvature; see for 
instance IIKKS89I |Kus9 1 1 IGKS03I . The rotational analog is the following torque balance; 



(p XX ijds — (p X X {ly X dx) —II H{x x ly) dA = / / xxH dA. 
J-Y Jy JJd JJd 

hat related is the following equation; 

cp X- rjds ^ cp X- {1/ X dx) = / / (H • x - 2) dA. 

J y J y JJd 



/ 7 J J JJd 

It gives, for example, an interesting expression for the area of a minimal (H = 0) surface. 



4. Discrete surfaces 

For us, a discrete or polyhedral surface AI C E'^ will mean a triangulated surface with 
a continuous map into space, linear on each triangle. In more detail, we start with an 
abstract combinatorial triangulation — a simplicial complex — representing a 2-manifold 
with boundary. We then pick positions p G E'^ for all the vertices, which uniquely deter- 
mine a linear map on each triangle; these maps fit together to form the polyhedral surface. 

The union of all triangles containing a vertex p is called Star(p), the star of p. 
Similarly, the union of the two triangles containing an edge e is Star(e). 
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4.1. Gauss curvature 

It is well known how the notion of Gauss curvature extends to such discrete surfaces M. 
(Banchoff fBa n67[ |Ban70ll was probably the first to discuss this in detail, though he 
notes that Hilbert and Cohn-Vossen IIHCV32I §29] had already used a polyhedral ana- 
log to motivate the intrinsic nature of Gauss curvature.) Any two adjacent triangles (or, 
more generally, any simply connected region in M not including any vertices) can be 
flattened — developed isometrically into the plane. Thus the Gauss curvature is supported 
on the vertices p G M. In fact, to keep the Gauss/Bonnet theorem true, we must take 



Here, the angles 6i are the interior angles at p of the triangles meeting there, and Kp is 
often known as the angle defect at p. If D is any neighborhood of p contained in Star(p), 
then §Qjj Kgds = ^Of, when the triangles are acute, this is most easily seen by let- 
ting dD be the path connecting their circumcenters and crossing each edge perpendicu- 
larly. Analogous to our intrinsic characterization of Gauss curvature in the smooth case, 
note that the circumference of a small e ball around p here is exactly 27r£ — eKp. 

This version of discrete Gauss curvature is quite natural, and seems to be the correct 
analog when Gauss curvature is used intrinsically. But the Gauss curvature of a smooth 
surface in E'' also has extrinsic meaning; for instance the total absolute Gauss curvature 
is proportional to the average number of critical points of different height functions. For 
such considerations, Brehm and Kiihnel IIBK82I suggest the following: when a vertex p 
is extreme on the convex hull of its star, but the star itself is not a convex cone, then 
we should think of p as having both positive and negative curvature. We let /\ + be the 
curvature at p of the convex hull, and set /v + — Kp > 0. Then the absolute 

curvature at p is /\ + + Kp, which is greater than \Kp\. This discretization, is of course, 
also based on preserving an integral curvature relation — a different one. (See also IIBK971 



We can use the same principle as before — preserving the Gauss/Bonnet theorem — 
to define Gauss curvature, as a measure, for much more general surfaces. For instance, on 
a piecewise smooth surface, we have ordinary K within each face, a point mass (again the 
angle defect) at each vertex, and a linear density along each edge, equal to the difference in 
the geodesic curvatures of that edge within its two incident faces. Indeed, clothes are often 
designed from pieces of (intrinsically flat) cloth, joined so that each vertex is intrinsically 
flat and thus all the curvature is along the edges; corners would be unsightly in clothes. 

Returning to polyhedral surfaces, we note that Kp is clearly an intrinsic notion (as it 
should be) depending only on the angles of each triangle and not on the precise embedding 
into E'^. Sometimes it is useful to have a notion of combinatorial curvature, independent 
of all geometric information. Given just a combinatorial triangulation, we can pretend 
that each triangle is equilateral with angles 9 = 60°. (Such a euclidean metric with cone 
points at certain vertices exists on the abstract surface, independent of whether or not it 
could be embedded in space. See the survey IITroOTII .) 




KdA := 2^ Kj 




lvDA95l .) 
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The curvature of this metric, Kp = -1(6 — deg\p), can be called the combinatorial 
(Gauss) curvature of the triangulation. (See IIThu981 llK"'"08l for combinatorial applica- 
tions of this notion.) In this context, the global form ^ Kp = 27rx(M) of Gauss/Bonnet 
amounts to nothing more than Euler's formula x = y — E + F . (We note that Forman has 
proposed a combinatorial Ricci curvature ||For03ll ; although for smooth surfaces, Ricci 
curvature is Gauss curvature, for discrete surfaces Forman's combinatorial curvature does 
not agree with ours, so he fails to recover the Gauss/Bonnet theorem.) 

Our discrete Gauss curvature Kp is of course an integrated quantity. Sometimes it is 
desirable to have instead a curvature density, dividing Kp by the surface area associated 
to the vertex p. One natural choice is Ap ;= -i Arca(Star(p)), but this does not always 
behave nicely for irregular triangulations. One problem is that, while Kp is intrinsic, de- 
pending only on the cone metric of the surface, Ap depends also on the choice of which 
pairs of cone points are connected by triangle edges. One fully intrinsic notion of the area 
associated to p would be the area of its intrinsic Voronoi cell in the sense of Bobenko and 
Springborn IIBS05I : perhaps this would be the best choice for computing Gauss curvature 
density. 

4.2. Vector area 

The vector area formula 

:= - / X X dx = ( ( vdA 

needs no special interpretation for discrete surfaces: both sides of the equation make sense 
directly, since the surface normal v is well-defined almost everywhere. However, it is 
worth interpreting this formula for the case when D is the star of a vertex p. More gener- 
ally, suppose 7 is any closed curve (smooth or polygonal), and D is the cone from p to 7 
(the union of all line segments pq for (767). Fixing 7 and letting p vary, we find that the 
volume enclosed by this cone is an affine linear function of p e E'^, and thus 

A If 

A„ := V„ Vol £> = — = - * X X dx 
3 6 

is independent of the position of p. We also note that any such cone D is intrinsically flat 
except at the cone point p, and that 2tt — Kp is the cone angle at p. 

4.3. Mean curvature 

The mean curvature of a discrete surface M is supported along the edges. If e is an edge, 

and e C D C Star(e) = Ti U T2, then we set 

He / / H dA = (p rj ds ^ e X — e X 1/2 ^ Jie — J2e. 

J Jd JdD 

Here vi is the normal vector to the triangle Ti, and the operator Ji rotates by 90° in the 
plane of that triangle. Note that |He| = 2 sin(0e/2) |e|, where 6e is the exterior dihedral 
angle along the edge, defined by cos^g ~ vi ■ V2. 
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No nonplanar discrete surface has Hg = along every edge. But this discrete mean 
curvature can cancel out around vertices. We set 



The area of the discrete surface is a function of the vertex positions; if we vary only one 
vertex p, we find that Vp Area(M) = —Hp. This mirrors the variational characteriza- 
tion of mean curvature for smooth surfaces, and we see that a natural notion of discrete 
minimal surfaces is to require Hp = for all vertices ||Bra92[|PP93| . 

Suppose that the vertices adjacent to p, in cychc order, are pi, . . . , p„. Then we can 
express Ap and Hp explicitly in terms of these neighbors. We get 



where and Pi are the angles opposite edge ppi in the two incident triangles. This latter 
equation is the famous "cotangent formula" IIPP93I IWarOSI which also arises naturally in 
a finite-element discretization of the Laplacian. 

Suppose we change the combinatorics of a discrete surface M by introducing a new 
vertex p along an existing edge e and subdividing the two incident triangles. Then Hp in 
the new surface equals the original Hg, independent of where along e we place p. This 
allows a variational interpretation of Hg. 

4.4. Minkowski mixed volumes 

A somewhat different interpretation of mean curvature for convex polyhedra is found in 
the context of Minkowski's theory of mixed volumes. (In this simple form, it dates back 
well before Minkowski, to Steiner MSte40l .') If X is a smooth convex body in E'^ and 
Bt{X) denotes its t-neighborhood, then its Steiner polynomial is: 



Here, by Gauss/Bonnet, the last integral is always An. 

When X is instead a convex polyhedron, we already understand how to interpret 
each term except J^^. H dA. The correct replacement for this term, as Steiner discovered, 
isJ2e^e |e|. This suggests iJg := 6*6 |e| as a notion of total mean curvature for the edge e. 

Note the difference between this formula and our earlier |He| = 2s'm{9e/2) \e\. 
Either one can be derived by replacing the edge e with a sector of a cylinder of length |e| 
and arbitrary (small) radius r. We have // H dA = He, so that 





and similarly 







HdA |He| = 2sin(6le/2) |e| < ee\e\^H, 




HdA. 
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The difference is explained by the fact that one formula integrates the scalar mean curva- 
ture while the other integrates the vector mean curvature. Again, these two discretizations 
both arise through preservation of (different) integral relations for mean curvature. 

See ISulOS l for a more extensive discussion of the analogous situation for curves: 
although as we have mentioned, the sum of the turning angles 'ipi is often the best notion 
of total curvature for a polygon, in certain situations the "right" discretization is instead 
the sum of 2 sin -0^/2 or 2 tan ipi/2. 

The interpretation of curvatures in terms of the mixed volumes or Steiner polyno- 
mial actually works for arbitrary convex surfaces. (Compare llSchOSl in this volume.) 
Using geometric measure theory — and a generalized normal bundle called the normal 
cycle — one can extend both Gauss and mean curvature in a similar way to quite general 
surfaces. See IFed59l lFu94l ICSM03ilCSM06l . 

4.5. Constant mean curvature and WUlmore surfaces 

A smooth surface which minimizes area under a volume constraint has constant mean 
curvature; the constant H can be understood as the Lagrange multiplier for the constrained 
minimization problem. A discrete surface which minimizes area among surfaces of fixed 
combinatorial type and fixed volume will have constant discrete mean curvature H in the 
sense that at every vertex. Hp = HAp, or equivalently Vp Area = — i/VpVol. 

In general, of course, the vectors Hp and Ap are not even parallel: they give two 
competing notions of a normal vector to the discrete surface at the vertex p. Still, 

_ I Vp Area I _ |Hp| _ | //st^..^^) H 
|VpVol| - |Ap| - |//s,,,(p)-dA| 

gives a better notion of mean curvature density near p than, say, the smaller quantity 

|Hp| _ iJJUdAl 
Ap JJldA ' 

where Ap := i Area(Star(p)). 

Suppose we want to discretize the elastic bending energy for surfaces, / j dA, 
known as the Willmore energy. The discussion above shows why J2p ^p^p (which was 
used in llFS+971 ) is a better discretization than J^p |Hp|^/Ap (used fifteen years ago 
in IIHKS92II ). Several related discretizations are by now built into Brakke's Evolver; see 
the discusion in IBraOTI . Recently, Bobenko IIBob05llBob08l has described a completely 
different approach to discretizing the Willmore energy, which respects the Mobius invari- 
ance of the smooth energy. 

4.6. Relation to discrete harmonic maps 

As mentioned above, we can define a discrete minimal surface to be a polyhedral surface 
with Hp = 0. An early impetus to the field of discrete differential geometry was the real- 
ization (starting with IIPP93II ) that discrete minimal surfaces are not only critical points for 
area (fixing the combinatorics), but also have other properties similar to those of smooth 
minimal surfaces. 
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For instance, in a conformal parameterization of a smooth minimal surface, the co- 
ordinate functions are harmonic. To interpret this for discrete surfaces, we are led to the 
question of when a discrete map should be considered conformal. In general this is still 
open. (Interesting suggestions come from the theory of circle packings, and this is an area 
of active research. See for instance MSteOSI IBH031 IBobOSI lKgS06l |Spr06| . ) 

However, we should certainly agree that the identity map is conformal. A polyhedral 
surface M comes with an embedding Id a/ : il/ ^ E'^ which we consider as the identity 
map. Indeed, we then find (following IIPP93I ) that M is discrete minimal if and only 
if Ma/ is discrete harmonic. Here a polyhedral map / : M ^ E'^ is called discrete 
harmonic if it is a critical point for the Dirichlet energy, written as the following sum over 
the triangles T of M: 

E{f) :=^|V/TpArcaM(r). 

T 

We can view E{f )~ Area f{M) as a measure of nonconformality. For the identity map, 
E{ld]\f) = Arca(Af) and Vpi?(Idj\f ) = Vp Arca(A/), confirming that M is minimal if 
and only if Id a/ is harmonic. 



5. Vector bundles on polyhedral manifolds 

We now give a general definition of vector bundles and connections on polyhedral mani- 
folds; this leads to another interpretation of the Gauss curvature for a polyhedral surface. 

A polyhedral n-manifold P" means a CW-complex which is homeomorphic to an 
n-dimensional manifold, and which is regular and satisfies the intersection condition. 
(Compare (ZieOS) in this volume.) That is, each n-cell (called a. facet) is embedded in P" 
with no identifications on its boundary, and the intersection of any two cells (of any di- 
mension) is a single cell (if nonempty). Because P" is topologically a manifold, each 
(n — l)-cell (called a ridge) is contained in exactly two facets. 

Definition. A discrete vector bundle of rank k over P" consists of a vectorspace 
Vf = E*^ for each facet / of P. A connection on is a choice of isomorphism 0,, 
between Vf and Vf for each ridge r = / n /' of P. We are most interested in the case 
where the vectorspaces Vf have inner products, and the isomorphisms <j)r are orthogonal. 

Consider first the case n = 1, where P is a polygonal curve. On an arc (an open 
curve) any vector bundle is trivial. On a loop (a closed curve), a vector bundle of rank k 
is determined (up to isomorphism) simply by its holonomy around the loop, an automor- 
phism (f>:E'' ^ E''. 

Now suppose P" is linearly immersed in E'' for some d. That is, each fc-face of P is 
mapped homeomorphically to a convex poly tope in an affine fc-plane in E'', and the star 
of each vertex is embedded. Then it is clear how to define the discrete tangent bundle T 
(of rank n) and normal bundle N (of rank d — n). Namely, each Tf is the n-plane parallel 
to the affine hull of the facet /, and Nf is the orthogonal {d — n)-plane. These inherit 
inner products from the euclidean structure of E''. 
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There are also natural analogs of the Levi-Civita connections on these bundles. 
Namely, for each ridge r, let e [0, tt) be the exterior dihedral angle bewteen the 
facets fi meeting along r. Then let 4>r '■ E'^ be the simple rotation by this angle, 

fixing the affine hull of r (and the space orthogonal to the affine hull of the fi). We see 
that (f>r restricts to give maps Tf Tf and Nf Nf ; these form the connections we 
want. (Note that T © = E'' is a trivial vector bundle over P", but the maps 0^ give a 
nontrivial connection on it.) 

Consider again the example of a closed polygonal curve C E'^. The tangent 
bundle has rank 1 and trivial holonomy. The holonomy of the normal bundle is some 
rotation of E''^^. For d = 3 this rotation of the plane is specified by an angle equal 
(modulo 2tt) to the writhe of P^. (To define the writhe of a curve as a real number, rather 
than just modulo 27r, requires a bit more care, and requires the curve P to be embedded.) 

Next consider a two-dimensional polyhedral surface P^ and its tangent bundle. 
Around a vertex p we can compose the cycle of isomorphisms (pe across the edges in- 
cident to p. This gives a self-map (pp -.Tf ~> Tf. This is a rotation of the tangent plane by 
an angle which — it is easy to check — equals the discrete Gauss curvature Kp. 

Now consider the general case of the tangent bundle to a polyhedral manifold P" C 
E''. Suppose p is a codimension-two face of P. Then composing the ring of isomorphisms 
across the ridges incident to p gives an automorphism of E" which is the local holonomy, 
or curvature of the Levi-Civita connection around p. We see that this is a rotation fixing 
the affine hull of p. To define this curvature (which can be interpreted as a sectional 
curvature in the two-plane normal to p) as a real number and not just modulo 2tt, we 
should look again at the angle defect around p, which is 27r — ^ (3i where the (3i are the 
interior dihedral angles along p of the facets / incident to p. 

In the case of a hypersurface P^^^^ in E'', the one-dimensional normal bundle is 
locally trivial: there is no curvature or local holonomy around any p. Globally, the normal 
bundle is of course trivial exactly when P is orientable. 
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